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We study exactly soluble system of trapped bosonic particles interacting by a model harmonic 
forces. The model allows for detailed examination of the order parameter (condensate wave function) 
as well as concept of the off-diagonal and diagonal order. We analyze the effect of interactions 
' on the condensate and show that sufficiently strong interactions, attractive or repulsive, lead to 

I destruction of the condensate. In the thermodynamic limit this destruction has a critical character. 

, ^ I It is shown that the existence of the coherent state of bosons is related to existence of two length 

^ ' scales determined by one- and two-particle reduced density matrices. The condensate can exist only 

if the two length scales are of the same order. Interactions, both repulsive and attractive, change 
their relative size which may lead to destruction of coherence in the system and depletion of the 
' condensate. We suggest that this scenario is model independent. 
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I. INTRODUCTION 



Recent advances in the trapping and cooling techniques which led to the achievement of the Bose-Einstein conden- 
sation of dilute gases have renewed interest in various aspects of many body theory. A cloud of weakly interacting 
trapped atoms is an ideal system for which various aspects of many body theory can be tested and verified. The 
Q ideal bosonic gas undergoes the Bose-Einstein condensation if the phase-space density exceeds one. This phenomenon 
manifests itself by the macroscopic occupation of the single particle ground state. In the case of an interacting system 
the condensate wave function can be defined by spectral decomposition of the one-body reduced density matrix. This 
I ' decomposition is closely related to the off-diagonal long range order and existence of the order parameter i.e. 

^ . the 'classical' field with given amplitude and phase commonly used in the theory of superfluidity |^,^ . Realization 
OO I of decomposition procedure is practically impossible because it requires a full solution of the many-body problem. 
Instead mean-field approaches are commonly used. The basic idea for a mean-field description of the dilute, weakly 
interacting Bose gas below transition temperature was introduced by Bogoliubov [Q . Most of the results for the inter- 
, acting Bose-Einstein condensate are obtained within the Bogoliubov theory which in many cases provides a reliable 
' quantitative description of the quantum Bose gas. On the contrary, the rigorous treatment is possible only for model 
: interactions. 

There are only few exactly soluble models of quantum systems where the interactions between atoms is chosen in 
the form allowing for the exact analytic solution. These are: (i) the one-dimensional model of impenetrable bosons 
introduced by Girardeau ||], (ii) its contact potential version formulated by Lieb ||]; (iii) the model of particles 
. interacting by harmonic forces [p^|-p^. Although in the first two cases the formal solution is given but in practice 
(~| ' the problem is still quite complicated and quantitative calculations can be done for a very small number of particles 
O ■ only The latter case seems to be much simpler because, as it has been shown in |l4|, it can be reduced to the 
^ ' problem of noninteracting particles in a harmonic trap. Therefore in the following we are going to examine, within 
J> this exactly soluble model, various concepts related to the interacting Bose-Einstein condensate. 

This paper is organized as follows. In Sec. II we present exact results for harmonically interacting bosons trapped 
within harmonic potential |Q. In Sec. Ill we find the analytic expression for the order parameter and study the effect 
■ of quantum depletion of the condensate as well as quantum fluctuations at zero temperature. In Sec. Ill we analyse 
the off-diagonal and diagonal order and show how off-diagonal long range order disappears when interactions become 
very strong. We finish in Sec. IV with some concluding remarks. 



II. DIAGONALIZATION OF THE HAMILTONIAN 

In our previous paper [ p^ we have shown the algebraic method of diagonalization of the Hamiltonian describing a 
system of many particles interacting via harmonic forces. The system under consideration consists of many particles 
confined by an external harmonic potential interacting by harmonic forces, i.e. two body interaction potential has 
the form: 
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v{^,-^,) = ^n\-K,-^,)\ (1) 

where defines the interaction strength and ct = +1 signifies the attractive interaction of particles placed at positions 
Xi and Xj whereas cr = — 1 - corresponds to repulsive interactions. The total Hamiltonian of the A'^-particle system 
has therefore the following form: 

^ 1 

i—l i<j 

Let us first recall some exact results of [Q. For the sake of simplicity we denote the set of all particle positions 
vectors by X^v — (xi, . . . , x^v)- The Hamiltonian which is a quadratic form of positions and momenta of particles can 
be easily diagonalized if one introduces collective variables: 



where = (xj, . . . ,x^) and the matrix Qjy = {Qij} is orthogonal. One of these collective variables namely the 



X^-Qa^X^, (3) 
center of mass of iV-particle system plays a particularly important role: 

AT 
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^E-^- (4) 
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The choice of — 1 remaining variables X^ -^ — (x^, . . . , x^_.^) is not unique but this does not lead to any physical 
implications. In particular: 



N-l N 



i=l 1=1 

In the following we are going to use a similar notation for description of a subsystem of s-particles, s — 1, . . . , N . 

The transformation defined above brings the Hamiltonian to the diagonal form and its eigenenergies can be easily 
found. However, while determining a spectrum one must also take into account the proper symmetry of a total wave 
function. In the case of bosonic particles {N > 2) the allowed energies are: 

E=('^+m) + C^{N-l)+n)u;, (6) 
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where m — 0,1,2..., n — 0,2,3... and w — Vl + crNil^. The first term describes excitations of the center of 
mass - d-dimensional harmonic oscillator of frequency equal to one. The second term in the Eq.(||) corresponds to 
excitations of A'^ — 1 relative degrees of freedom. The frequency uj characterizes the effective potential experienced by 
— 1 collective modes. 

Let us observe that lu = I corresponds to the noninteracting case, the repulsive interactions give < lu < 1 while 
attractive forces lead to > 1. Moreover, very small values of w « signify very strong repulsion which almost 
destabilizes the whole system. It is very convenient to parameterize to by an exponent k defined in the following way: 

UJ = N". (7) 

This exponent can be related to the actual strength of the interaction. In fact, for weakly interacting gas (a; « 1) 
we obtain very small values of this parameter: k « 0, while for strong interactions « - repulsion, a; ^ 1 - 
attraction) we have |k| 3> 1. Moreover, k is positive in the case of attraction while it is negative for repulsion. Let us 
add at this point that in realistic situations of short-range interparticle interactions, large Bose-Einstein condensates 
can exist only for repulsive forces. In the case of attraction the size of the trapped condensate is limited to about 
1500 atoms In our oscillatory model the forces between particles are negligible at small distances, therefore the 
model leads to the condensation (in the thermodynamic limit) in both attractive and repulsive case. 
The ground state of the system is the following: 

Vl/(XAr) = <i>o(V^X$v_i)$o(x^), (8) 
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where (X5i^_-^, x^^) = QatX^v and the function <I>o(\/^^Ar-i) corresponds to the ground state of a system of iV — 1 
independent quasi-particles (in d spatial dimensions) interacting with an external potential of the harmonic oscillator 
of frequency w: 

$o(%AJX$v_i) = (-) exp [-c.(X$^_i)V2] , (9) 

and $o(x^) is the ground state of the single particle (center of mass) trapped into harmonic potential: 

$o(x5v)= (^)'^'exp[-(x^)V2]. (10) 

Construction of excited eigenstates is difficult because it is not easy to impose the desired symmetry on the wave 
function. This procedure was describe in details in p4[|. 



III. ORDER PARAMETER AND QUANTUM DEPLETION 

If the energy of the system (or equivalently the temperature) is sufficiently small we expect that the system forms a 
Bose-Einstein condensate. The BEC of the ideal gas manifests itself by a macroscopic occupation of the single particle 
ground state. In the case of interacting system it is not obvious what is this particular state which is 'macroscopically 
occupied'. The identification of the macroscopically occupied quantum state is equivalent to the definition of the order 
parameter - the single particle wave function which is inherently related to the Bose condensation. The condensate 
subsystem can be then quite accurately described by the A'o-fold product of the order parameter, where Nq ~ 0{N) 
is the occupation of the condensate. 

In the following we use our model to demonstrate how to define the order parameter, occupation of the condensate, 
and its fluctuations. At zero temperature the system is in the ground state and one might naively expect that it 
is totally Bose condensated. However, the ground state of the N-particle bosonic system is not equivalent to the 
Bose-Einstein condensate. Interactions can significantly deplete the condensate. We are going to show this effect in 
the most spectacular but also in relatively simple case of the zero temperature. 

Let us now define the hierarchy of the reduced s-particle density matrices which can be conventionally obtained by 
averaging the density matrix of the total system of N particles over the degrees of freedom oi N—s remaining particles. 
For a given iV-particle quantum state ^(X.n) the corresponding s-particle reduced density matrix ps(X.s; Y^) is defined 
by: 

Ps{X,;Ys)= y"dRjv-s**(X„RAr_,)*(Y„RAr_,). (11) 

We use previously defined shorthand notation for vectors in the configuration space of s-particles, for example Xjv = 
(Xs, Rjv_s) and X^ = (xi, . . . ,Xs). The reduced density matrix describes the subsystem of s-particles and can be 
directly related to different measurement processes. For the statistical description of the system one should first of 
all define the statistical density matrix by averaging all TV-particle density matrices with the appropriate statistical 
weights depending on the ensemble. In general it is quite a complicated task but at zero temperature there is only 
one quantum state of the system and no statistical averaging is necessary. 

The total wave function (or density matrix) carries all the information about the system. In real experiments 
however one does not probe simultaneously all the particles. Typical detection scheme consists on the measurement 
of one or at most few particles at a given time. In other words a single measurement process is reduced to a subsystem 
of small number of particles. Such subsystems are described by reduced density matrices. In the considered here case 
of zero temperature the s-particle density matrix can be brought to the following form: 

p,(X,; Y,) = p^^^(x^,y,^)<^>o(^A3X^_l)<i>o(V:^Y:_l)■ (12) 
The functions $o(v^Xg_;^) describes the ground state of s — 1 quasi-particles: 

d(s-l)/4 

MV^^U) = (-) exp [-c.(X^_i)V2] , (13) 

while p'^^^(x^,y^) corresponds to the density matrix of center of mass of the subsystem: 
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The s-particles collective coordinates are defined in the familiar way: (X3_]^,X3) = Q^Xs, (Y^_]^,y^) = Qs Ys and 
frequencies ajg, Sg as well as auxiliary parameter 7^ are: 

(15) 



TV 



7s 



1 -ojV s{N 



(16) 

Having defined the s-particle matrices we are ready now to analyze the nature of the Bose-Einstein condensation 
of the interacting system and to discuss the meaning of the order parameter. To this end we write the density matrix 
Eq.(|l2|) in the diagonal form: 

p,(X.;Y,) = ^A(f)(/)(^)(X,)0(f)(Y,). (18) 

n 

The function (/'[^•'(Xs) can be treated as the wave function of the s-particle subsystem: 

0(fHx,) = $o(V^X^_i)$n(V^x^), (19) 

where $o(v^Xs-i) is the ground state wave function of the relative degrees of freedom. This function corresponds 
to the ground state of s — 1 noninteracting quasi-particles (in d-spatial dimensions) subject to the external harmonic 



potential of frequency ui. The second part of the Eq.(19) describes states of the center of mass of s-particles; is 
simply the d-dimensional harmonic oscillator wave function corresponding to the effective center of mass frequency 
as ■ Quantum numbers n = (m , . . . , n^) label different states of the center of mass while n = ni + . . . + nd corresponds 
to the energy of the given state. The effective center of mass frequency as is: 

as = [iUsius + 6s)]^^\ (20) 

It is interesting to observe that all the frequencies of the relative motion of the s-particles subsystem are exactly the 
same as the frequencies of the relative motion of the whole system, i.e. equal to ui. On the other hand the center 
of mass oscillation frequency of the subsystem is neither equal to co nor to 1 (trap frequency). This collective degree 
of freedom couples to the center of mass of iV — s remaining particles what leads to some energy shift. Finally, the 

(s) 

eigenvalues A„ of ps are equal to the occupation probabilities of a given s-particle state: 

- f^sV^' pyZ^ y / as-Us V _ ^21) 
as J \ujs + asj \as+u!sj 

It follows from the normalization condition for the density matrix that: Xn^ ~ 1. 

The spectral decomposition of the reduced one-particle density matrix gives natural single-particle states 0n^(x). 
These states are crucial for the definition of the condensate wave function (order parameter). It can be seen from 
Eq.(|l]) that if N goes to infinity (thermodynamic limit) with fixed value of the interaction strength k the lowest 
eigenvalue Aq^-* dominates the others and in the limit of weak interactions we get: 

A« ^ 1, (22) 

A«-(^^)", if n^O. (23) 

This behavior signifies nothing else but the onset of the Bose-Einstein condensation. The single particle density 
matrix becomes very close to the pure state because with quite good accuracy it can be approximated by pi (x, y) « 
(/)o^^(x) 0Q^^(y). This particular single-particle ground state (/)g^^(x) is usually called the order parameter. The A^- 
particle wave function can be quite accurately approximated by the A^-fold product of the order parameter. 
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Our analytic formula allows to study quantitatively the role of interactions on the Bose-Einstein condensate. On 
the basis of the discussion it is obvious that the average occupation of the condensate becomes: 

{No) =N j dxdy 0« (x) pi (x, y) (y) - N\^^^ . (24) 

In the case of the ideal gas at zero temperature the above equation gives, of course, (A^o) — all particles occupy 
the single particle ground state. For studied here model of the interacting system it is more convenient to use as 
independent variables the pair of k and N rather than lu and N. The exponent k occurs to be inherently related to the 
oscillatory interactions because when using this parameter we discover some universal features. For a fixed number 
of particles, if the interaction strength |k| ~ logo; grows, the occupation of the condensate decreases. This behavior 
is presented in Fig. 1 where we show the mean occupation of the condensate versus the exponent k = logo;/ log 
for different values of particle number N in three spatial dimensions {d = 3) . The values of k less than zero signify 
repulsive interactions while k > corresponds to attraction. One can easily see that if the interaction becomes strong 
(|«;| ~ 1) the condensate is almost totally depleted. All curves presented in the figure tend to an universal curve if 
the number of particles increases. When A^ increases to infinity with k being constant then our expression for the 
occupation of the condensate has the form: 

(^o) _( 2 
N Vl + VN''-^ + N'^^+^^ + 1 

The above formula, valid in the thermodynamic limit, gives an universal critical behavior. It exhibits no depletion 
(A'o = N) for \k\ < 1 followed by an abrupt jump and total destruction of the condensate (A'o = 0) for |k| > 1. 
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FIG. 1. Mean occupation of the condensate plotted as a function of the parameter k = log a;/ log A'' for different number of 
particles; A^ = 10^ - dashed line; A^ = lO'' - dotted line, and A^ = 10** - full line. 

The effect of quantum depletion of the trapped atomic condensate with a short range interactions, for the realistic 
experimental parameters, has been estimated to be of the order of 1% Q. This is opposite to the case of superfluid 
helium where this effect accounts for depletion as big as more than 90% |l5|. Our model exhibits very interesting 
feature. It shows that in large A^ limit the quantum effects are almost negligible or totally destroy the condensate 
depending on the value of the interaction strength. At this point it is not clear if this is an unique feature of our 
model or if it is more general result. 

The 2-particle reduced density matrix allows to find a joint probability of finding a particle in a given single particle 
state provided that another particle is also in some given state. In particular we have: 



{Nq{No - 1)> = A^(A^ - 1) / dXadYa 



P2(X2, Y2) 



(y2)4'^(yi) 



(26) 



Simple integration gives: 



{No{No-l))=N{N-l) 



\ W2 + ai 



W2 + 0L\ 

W2 + ai + ^2 



dj2 



(27) 
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Now we are ready to analyze the fluctuations of the condensate defined as: 

(S'No) = {{No)') - {Nof. (28) 

These fluctuations are shown in Fig. 2. We see that as the interaction strength grows (at fixed number of particles) the 
fluctuations start to grow from zero value for the ideal gas. However, again when the interactions become so strong 
that condensate practically disappears (|k| ~ 1) fluctuations also decrease - as there is no condensate the fluctuations 
also die out. The fluctuations are maximal in a region of the critical destruction of the condensate by quantum effects. 
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FIG. 2. Fluctuations of the condensate plotted as a function of the parameter k — logo;/ log for different number of 
particles; = 10^ - dashed line; = lO'' - dotted line, and N = 10^ - full line. 



IV. OFF-DIAGONAL ORDER AND CORRELATIONS 



In this subsection we are going to study (within our model) the phase and density correlations for the ground state 
of the interacting system of bosonic particles. In the early works of Penrose and Onsager [Q it has been shown that 
phenomenon of Bose condensation is evidenced by the presence of the off-diagonal long range order in the one-particle 
density matrix pi (x, y) : 

lim pi(x,y)^0. (29) 

x-y|^oo 

The above condition signifies a large scale correlations of the off-diagonal elements of the single-particle matrix. The 
limiting value of correlations is called the off-diagonal long range order parameter. In a finite system this concept 
cannot be strictly applied as the reduced density matrix vanishes at large distances. However, for a practical purposes 
the off-diagonal long range order can be related to the behavior of the off-diagonal elements of the density matrix at 
distances equal to the size of the system. 

In order to identify long range correlations we should first analyze characteristic length scales. It is natural to 
relate the size of the system to the spatial extension of the single particle density, i.e. the diagonal elements of the 
one-particle matrix: 

ni(x) = /^pi(x,x) = N (^'^y^^ exp{-ujix'). (30) 



Therefore, the size of our system is: 



(31) 



For large values of N an extension of the system for attractive interactions approaches the value of one, £i w 1, i.e. 
is governed by the characteristic length scale of the external potential felt by the center of mass. Relative degrees of 
freedom are localized within much smaller distances 1 / ^/lu <^ £i . For repulsive interactions a spatial extension of the 
system is related to the length scale of a mean field potential experienced by relative degrees of freedom, £i « Xj ^fuj. 
In this case the magnitude of the center of mass spreading is relatively small 1 <C £i. These features have serious 
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implications on behavior of the order parameter and ofF-diagonal order. Let us stress that existence of different length 
scales - the one related to the external potential, and another related to the interparticle interactions is the general 
feature of all two-body central forces. This is not an unique behavior of considered here long range oscillatory model. 

The presence of the ofF-diagonal long range order in the case of inhomogeneous system of finite spatial extension 
manifests itself by a large distance behavior of the following correlation function: 



5i(x,y) 



y/ni{x)ni{y) 

Accordingly, we shall refer to notion of the off-diagonal long range order if: 

lim gi(x,y) ^ 0. 
|x-y|->^i 



(32) 



(33) 



By dividing the correlation function by the square root of single particle densities we ensure that correlations are 
not sensitive to the local density of particles. Let us observe that off-diagonal correlation function describes phase 
correlations in our system. In the case when the one-particle density matrix corresponds to a pure state ^(x) = 
i?(x) expos' (x)] (where i?(x) and S{x) are the modulus and phase of the wave function), then (?i(x, y) — exp{— i[S'(x) — 
5'(y)]}. In the context of the Bose-Einstein condensation of trapped alkali gases such phase correlations in a trapped 
2-dimensional interacting system has been studied recently in Ref. |16|. 
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FIG. 3. Off-diagonal correlations gi{x,0)) for a system of A*' = 10^ particles for different interaction strength 
K = log a;/ log TV. Note that gi{x,0) does not depend on the sign of interaction. The full lines correspond to the follow- 
ing values of |k| (from the top to the bottom); |k| = 0.3, 0.6, 1.3. The dashed line corresponds to the single particle density 
function n\(x)/N. The distance has been scaled in such a way that single particle density does not depend on the strength of 
interaction. 



Now we are going to apply the above definition to study the coherence of the ground state of bosonic particles 
interacting by harmonic forces. Similarly as in the previous section we will investigate the large TV limit varying the 
interaction strength k. Because the system is inhomogeneous the correlation function depends on the six variables 
(in three dimensional space). To simplify our analysis we fix the position y at the trap center, y = 0. With this 
restriction the correlation function depends only on the distance x = |x|: 



5i(a:,0) = exp 

what gives for the limit of large distances, x ^ £i, the expression for the off-diagonal order: 

,,(f„0) = exp(-A 



(34) 



(35) 



In the Fig. 3 we show the behavior of the off-diagonal order Eq.(|34|) as a function of distance to the trap center for 
N — 10^ particles and different interaction strength k. We see that for a weak interaction the system exhibits the 
off-diagonal long range order, while if interactions are strong the phase correlations rapidly vanish at distances smaller 
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than the natural size of the system. For comparison we show in the Fig. 3 the single particle density n{x)/N (dashed 
line). This kind of dependence of the off-diagonal order is consistent with the behavior of the order parameter (the 
condensate wave function) discussed previously. Indeed, the coherence length L^^ff which appears in the equation 
Eq.(^5[) in the large N limit has a form: 

-(1) _ / CJi ^ 



^"ff~[sj - v^— 1 + ^^^^ 

Let us notice that the coherence length is governed by the same type of K-dependence as the mean occupation of the 
condensate Eq.(p5|). Therefore it is not surprising that the order parameter and off-diagonal order behaves similarly. 
Disappearing of the off-diagonal long range order has exactly the same critical character as destruction off the order 
parameter when the interaction strength exceeds value one, \k\ > 1. The coherence length in large N limit, N ^ oo 
is: 

4// = oo, if |«;|<1, (37) 

4// = 0, if \i^\>l- (38) 

So far we have studied the off-diagonal correlations for the one-particle reduced matrix. However, a similar analysis 
can be done for any s-particle density matrix. Without going into unnecessary algebra one can easily estimate the 
coherence length for a given s-particle subsystem. The form of the corresponding reduced density matrix (see Eq.(p^) 
explicitly suggests that the order of magnitude of the off-diagonal elements at distances comparable to the size of the 
system, x = £i can be estimated by: exp[— 5s/(4ws)(a;/£i)^]. Therefore, the s-order coherence length is: 



UJs 



1/2 



" IZ ) " V^(7v--i + iv-i)- ^^^^ 

This expression shows that there is an off-diagonal long range coherence on each level of hierarchy of the reduced 
density matrices provided that |k| < 1. 

This statement, however, requires some comment. In fact the equation Eq.(|39|) is valid only if \imN^oois/N) = 0, 
i.e. if we consider the subsystem of finite number of particles (while the total number of particles is infinite). On the 
contrary, in the limit limN^oois/N) = (3 > we have: 



UJs 



1/2 



if ) = ( ^ ) = ^ 0. (40) 

There is no coherence in large subsystems when the number of particles in the subsystem is of the order of the total 
number of particles, s ~ 0{N). However, as long as the detection scheme is reduced to an observation of at most 
few particles simultaneously (what is true for all experiments) the Bose system at zero temperature can be viewed 
as described by the coherent state of many particles. This point of view is justified on every level of the reduced 
s-particle density matrix provided that s/N ^ 1. This kind of behavior is in a spirit of Yang ( ^) conjecture, who 
suggested that Bose-Einstein condensation is evidenced by the fact that every s-particle matrix can be related to a 
simple s-fold product of the order parameter. In our model of the interacting condensate the above statement is true 
for a whole hierarchy of s-particle matrices when s/N — > 0. If s/N <^ 1 and |k| < 1 simple integration gives: 



where /3'^*-' is: 



p.(X,, Y,) ^ ■ ■ ■ <\^s)] K'^(yi) . . . ^^o\ys)], (41) 



\LJ + ai J \ujs + aij \Us + ai + Ss J 

The equation Eq.(|4l|) is not valid if one considers s-particle matrices for s being of the order of the total number of 
particles, s ~ 0{N). At this point we want to make a comment about the notion of the coherence of the Bose-Einstein 
condensate. Approximate methods assume explicitly that a mean value of the boson field operator is different than 
zero in the case of the Bose-Einstein condensate. Therefore, the folk wisdom associates the condensate with the 
coherent state - the analogue of the coherent state of the electromagnetic field. This analogy is of limited value and 
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in fact may be misleading because the condensate must be in a Fock state in which a mean value the field operator 
vanishes. However, there is coherence in the condensate in the sense that majority of particles is described by the same 
wave function with the same phase. The coherence of the condensate manifests itself by the off-diagonal long range 
order. The coherence length of a phase of the condensate is at least of the same magnitude as the spatial extension 
of the system. Moreover this coherence is present also when one considers the subsystems of any number of particles 
significantly smaller than a total particle number. This facts allows for substantial simplification of description of 
many particle system. When the system is coherent, i.e. when one can assign exactly the same wave function to every 
particle then a behavior of the condensate can be correctly described by this single-particle function. The above does 
not apply to measurements based on simultaneous detection of almost all particles. In such a case the A''-body wave 
function has to be used instead. 

The presented analysis of the order parameter and off-diagonal correlations exposed a kind of symmetry with 
respect to change of the interaction sign. The differences can be seen while studying the diagonal correlations of 
the reduced density matrices. The diagonal elements of one-particle matrix are related to single particle density, 
ni(x) = 7Vpi(x, x) while p2(x, y;x, y) describes a joint probability of finding a particle at position x provided than 
another particle is at position y. The corresponding two-particle density (if x ^ y) is: 



n2(x, y) = N{N - l)p2(x, y; x, y). 
The above function can be related to two-particle correlation function: 

n2(x,y) - ni(x)ni(y) 



52(x,y) 



ni(x)ni(y) 



(43) 



(44) 



If two-particle density can be written as a product of corresponding one-particle ones then there are no diagonal 
(density) correlations, i.e. 5'2(x, y) = 0. Positive values of signify some clustering effect while negative values 
correspond to some effective repulsion. Note that g2 can be in principle infinitely large (extremely strong attraction) 
but its minimal value is limited by —1. This value corresponds to extremely strong repulsion when particles tend to 
avoid each other. 

In the case of studied here model, all density functions are inhomogeneous. Therefore, for the sake of simplicity we 
again fix one of the position vectors at the trap center, i.e. we set y = 0. Then the correlation function g2 depends 
only on the distance x = |x|: 



52 (a;, 0) = 



N 



1 / UJU12 \ 



d/2 



exp 



-- (w -h W2 - 2a;i) x 



1. 



The above equation gives the characteristic length scales of two-particle diagonal correlation: 

1/2 



r(2) 
diag 



U>1 



■ L02 — ^LVi 



(45) 



(46) 
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FIG. 4. Diagonal correlation g2{x,0)) for a system oi N = 10^ particles for different interaction strength k ■ 
K = 0.4 - dashed line; k = 0.8 - dotted line; k = 1.1 - full line. 
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Let us first consider the attractive forces. In this case the diagonal correlation function strongly depends on the 
magnitude of the interaction strength For weak interactions, < k < 1, the diagonal correlation length becomes 
very large L^"^^ — > oo similarly as the off-diagonal one. This is however, not necessarily the signature of strong 
correlations in the system. In the considered case the magnitude of correlations vanishes in large N limit within the 
whole system, 32(2;, 0) « if \x\ < £1. If attraction becomes very strong, k > 1, then the diagonal correlation length 
drastically decreases: 

41-l/V^=4?/-0- (47) 

In fact it becomes equal to the phase coherence length. Simultaneously the magnitude of these diagonal correlations 
grows to infinity, 52(2;, 0) « (jV""^)''/^ if |a;| < ii. This short range diagonal order is related to the clustering of 

particles. The joint probability of finding two particles very close to each other is very large provided that their 

(2) 

separation is smaller than L^^.^^g. Although the particles are clustered the position of this 'clustering spot' is localized 

(2) 

within the large distance determined by the spatial extension of the center of mass ground state, Lj^^^^g ^ £1 ~ 1. This 
behavior is illustrated in Fig. 4. The detailed comparison with corresponding off-diagonal correlations. Fig. 3 is very 
interesting. The presence off the off-diagonal long range order in the system is accompanied by the absence of the 
diagonal order for relatively weak interactions, k < 1. If interaction becomes very strong, k > 1, then a diagonal order 
becomes very large on a very small distance. Simultaneously the off-diagonal order disappears. Both characteristic 
length scales for diagonal and off-diagonal correlations are the same. 

For the repulsive forces the diagonal order depends very weakly on the interaction strength. The magnitude of 
diagonal correlations tends to zero with increasing particle number, g2{x, 0) « 2N~^ — SN*^ while the correlation 
length is large ij^^^g = iV/Vl ^ 2N'^ + 2N^^ 00. The change of the behavior of the diagonal order with the 
interaction strength k is similar to the case of attractive forces (see Fig. 3). However, as the correlations are practically 
negligible a small diagonal order which appears when k < — 1 is almost unvisible. In general we can say that similarly 
as in the case of attraction, there are no diagonal correlations when the system exhibits the off-diagonal long range 
order. 

V. CONCLUSIONS 

In our paper we used exactly soluble many-particle model to illustrate a rigorous procedure of dcfinig the condensate 
phase at zero temperature. By diagonalizing one-particle reduced density matrix we were able to study in details the 
role of interactions on the condensate. If the interaction strength becomes large |k| > 1 the condensate disappears 
even when the system is in its ground state. In the oscillatory model studied here this total depletion of the condensate 
has a critical character in the thermodynamic limit. 

The existence of nonvanishing order parameter is accompanied by the off-diagonal long range order. The off-diagonal 
order is present on all levels of the hierarchy of reduced s-particle density matrices {s/N <C 1) as long as |k| < 1. 
There is no diagonal long range order if the system exhibits the coherence (off-diagonal order). 

All these results have been obtained in the framework of a model. The model has some unrealistic features, its 
advantege is that it is exactly soluble. Therefore all the results presented in the paper are rigorous within the model. 

Most of the results are, however, model-independent and in fact were discussed previously but without rigorous 
proof. Thus we provided another way of justifying them by showing that they can be proved for model interactions. 

In particular the existence of (at least) two diffeent length scales, the one related to the external potential and the 
other resulting from the combined effect of interparticle interactions is a generic feature of all interacting systems. 
Moreover, as the first scale does not depend on the number of particles and interaction strength the second one does. 
The last scale is related to the spatial extention of the relative degrees of freedom. Our model shows that there is a 
coherence in the system if the two length scales are of the same order (different collective eigenmodes are not spatially 
separated). However, when these two scales become drastically differrent there is a spatial separation of different 
phases and off-diagonal long range order is destroyed. 

In the case of attraction the effect of destruction of the condensate relays on the fact that particles cluster in the 
region of the size much smaller than the center of mass range. On the contrary, for repulsive forces the condensate 
disapears when the amplitude of the center of mass zero-point oscillations is much smaller than the size of the system. 
In fact in the case of realistic attractive forces it is well known that interactions lead to clustering of atoms. Therefore 
condensates with negative scattering length are limited to relatively small number of atoms Our results suggest 
however, that there is a similar limitation on the size of the condensate interacting by repulsive forces. In fact in the 
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case of short range forces the Thomas-Fcrmi approach gives the ratio of the condensate size, R, to the characteristic 
length scale, aho, of the trap potential in the form: R/aho = {IBNas/ahoY^^ whre is the scattering length. For the 
largest existing trapped condensates {N ^ lO'') this value is about 10. This somehow fits perfectly to the scenario of 
destruction of the condensate by interactions presented in this paper. 
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